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Abstract The Cryogenic Dark Matter Search (CDMS) utilizes large mass, 3" di- 
ameter x 1" thick target masses as particle detectors. The target is instrumented 
with both phonon and ionization sensors, the later providing a~l V cm -1 electric 
field in the detector bulk. Cumulative radiation exposure which creates ^200 x 10 6 
electron-hole pairs is sufficient to produce a comparable reverse field in the detec- 
tor thereby degrading the ionization channel performance. To study this, the exist- 
ing CDMS detector Monte Carlo has been modified to allow for an event by event 
evolution of the bulk electric field, in three spatial dimensions. Our most resent 
results and interpretation are discussed. 
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1 Introduction 

The Cryogenic Dark Matter Search experimental US es a combination of phonon 
and ionization signals produced in radiation interactions in the detector target mass 
to determine event energy and interaction type. Electron recoils, which constitute a 
background for the experiment, result in 25% ionization energy with the remaining 
75% energy in the phonon system. Nuclear recoils, which make up the signal 
region, result in a reduced ionization signal to ~l/3 the amount for gamma events, 
with the remaining amount in the phonon system. Cuts are imposed on the ratio 
of ionization and phonon energy to allow passage of nuclear recoil events. 

While not the focus of this paper, events located outside of the bulk fiducial 
volume result in reduced ionization energy leading to electron recoils appearing 
like nuclear recoils. One example is surface events in which the high initial kinetic 
energy of charge carriers causes them to be injected into the wrong electrode. 
Another example are events at large radius, where fringing fields causes charge 
carriers to transport into the detector sidewalk 

There are also detector operating states which result in reduced ionization en- 
ergy readout, again causing misidentification of electron recoils as nuclear recoils. 
This state is described by a build up of space charge in the detector which sets up 
an electric field opposite to the applied field. This results in a net zero electric 
field which does not accelerate charges into the readout channel and causes an 
increased charge trapping rate within the detector bulk. 

2 Charge Transport 

Charge transport in germanium at low temperatures is described in more detail 
in another paper in these proceedings and elsewhere 34 . In our CDMS detector 
Monte Carlo, transport stops when either a charge traps on impurities or reaches 
the detector- vacuum boundary. While there is variation in trapping times between 
the CDMS detectors, the Monte Carlo models used in this paper used a represen- 
tative value deduced from calibration data, 8 /is and 64 jis for electron and hole 
trapping respectively. The carrier stopping location is then recorded and a pro- 
cess of updating the electric potential inside of the detector, to include the initial 
potential solution and the contribution from these charges begins. 

A component of the potential evolution calculation involves solving the elec- 
tric potential kernel of point charges inside the detector. Boundary conditions 
that describe the kernel include condition V = on the lithographically defined 
metal surfaces and conditions £ 0Ut dV 0Ut /dn — £i n dVi n /dn — —o and dVi n /dp — 
dV 0Ut /dp, where dn and dp represent differential lengths in the normal and paral- 
lel directions, at the germanium- vacuum boundary. Due to the complicated shape 
of the metal surfaces this is a non-trivial kernel to calculate and simpler approxi- 
mations are used for two different detector types. For the iZIP style detector^, the 
top and bottom of the detector have a relatively sparse 6.1% aluminum area cov- 
erage and are spaced ^2 mm away from the neighboring detector resulting in a 
topology that can be described as a long cylinder of germanium when vacuums are 
ignored (considered to have a Ge dielectric of £ — 16). There is a ground boundary 
condition provided by a copper supporting structure which surrounds the detec- 
tor. This structure is hexagonal and has a minimum center to copper distance of 



3 



38.7 mm, where the detector radius is 38.1 mm. With these and an additional ap- 
proximation that the copper can is described by a 38.7 mm radius cylinder, a soup 
can kernel is found to beJ^ 



V(p,0,z) = ^A nr 4(^rp)cos(^(0-0 o ))sinh(^ r (L + z)), z < zo (1) 

n=0r=0 

where 

A = 4(2 -8 n0 ) sinh(fc wr (L-z )) J n {k nr po) ^ 
a 2 sinh(2& nr L) k nr (J n j r \(k nr a)) 2 ' > 

and the Bessel function J n is defined to have zeros at the cylinder walls {r — a). The 
kernel below the source is found by interchanging + and — in the sinh function 
argument for both V and A. 

A high number of n terms constrains the solution in the angular direction and 
a high number of r terms constrains the solution in the radial and longitudinal 
directions. Taking the sums to infinity (or some large number) is computationally 
expensive given the 200k mesh points and 30 charge carrier pairs. Instead, the 
sums are computed to n = 3 and r = 81 for the iZIP detector resulting in a spread 
of ~ 71 / 2 in the angular direction and ^2 mm in the radial and longitudinal direc- 
tions. This spread is reasonable since the iZIP data was exposed to a high energy 
barium- 131 gamma source (356 keV) resulting in events distributed throughout 
the detector bulk. With the exception of A nr , all components in V are independent 
of charge location and for computational efficiency are computed in an initializa- 
tion procedure and stored in memory for reuse. Charge carriers in the next event 
are then propagated with this updated potential solution and the procedure is re- 
peated. 



3 iZIP Monte Carlo Results and Comparison to Data 

In addition to updating the electric potential, the charge carrier stopping position 
can also be run through an instrumentation simulator to recreate the detector re- 
sponse and allow a direct comparison between Monte Carlo and data. There are 
some measurement quantities which show qualitative similarities between Monte 
Carlo and data and some which do not show good agreement. 

An indicator of good agreement is seen in the hole radial-charge-partition 
quantities for the iZIP detector. This signal is defined to be the signal in the in- 
ner electrode minus the signal in the high radius electrode. As shown in Figure [T] 
both data and Monte Carlo show this measurement becoming negative as the field 
evolution progresses indicating a build up of holes at large radius near the hole 
collection electrode. 

A second indicator of good agreement is seen in the sum of electron and hole 
collection in the iZIP detector. As shown in Figure [2] both Monte Carlo and data 
have a reduced response at long times. A similar feature is seen in more detailed 
studies described in a companion paper in these proceedings. 7 This is generally 
assumed to be due to an induced electric field from trapped charges counteracting 
the applied electric field resulting in a reduced electric field magnitude in the de- 
tector. With this reduced electric field, the charges transport more slowly through 
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iZip Monte Carlo iZip Calibration Data 




2 4 6 8 10 12 14 o 0.2 o.4 o.6 o.s i 1.2 1.4 1.6 1.8 2 



Integrated Charge Carriers [number] x ^ 9 Integrated Charge Carriers [num, estimated] X 1 8 

Fig. 1 Radial charge partition from holes as measured in a CDMS iZIP detector. The figure 
on the left shows a Monte Carlo response of the detector whereas the one on the right is from 
calibration data. The vertical dashed line indicates periods of similar estimated exposure. 



the detector increasing their chance of becoming trapped and resulting in a re- 
duced charge collection signal. 

An indicator of poor agreement between Monte Carlo and data is seen in the 
charge partitioning along the longitudinal direction of the iZIP detector. The data 
in Figure [3] shows a gradual decline from to -0.1 indicating that electrons are 
transporting less efficiently to the charge electrode as the run progresses. This 
gradual decrease is strikingly absent in the Monte Carlo and instead there is a 
dramatic reduction of bulk events (partition = 0) and the measurement becomes 
increasingly noisy. 

A second indicator of poor agreement is observed in all of the measurements, 
notably the time scale for any change in measurement are about a factor of 100 
greater than what is seen in data. This suggests one of a few possibilities, that the 
detector was not initially in applied-field only condition when the dataset acquisi- 
tion was started, that some other mechanism other than gamma-ray interactions is 
generating charge carriers which reduce the electric field, or that the charge trap- 
ping rates are incorrect. In Figures [T] [2] and [3| the number of charge carriers in the 
Monte Carlo is known, however the number of carriers in data is inferred from the 
energy of events with energy >40 keV. It is additionally possible that events could 
have been created by lower energy, unresolvable events. 



4 Test Device Monte Carlo Results and Comparison to Data 

In addition to the iZIP detector, a charge transport test device with four concen- 
tric ionization channels patterned on one detector face was run. This test device 
is 100 mm in diameter and 33.3 mm thick and the top and bottom surfaces have 
a high ^100% area coverage. Additional discussion of this device and operating 
characteristics is found in these proceedings. 8 This detector is also surrounded by 
a hexagonal copper can with minimum center to can distance of 5 1 .5 mm. This de- 
tector was exposed to Am- 241 gammas (59.5 keV) from four collimated, surface 
sources. Since the sources (and therefore events) are more localized it is prefer- 
able to more tightly constrain the kernel in the angular direction and sums are 
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iZip Monte Carlo 




iZip Calibration Data 
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Fig. 2 Charge signal as measured in a CDMS iZIP detector. The figure on the left shows a 
Monte Carlo response of the detector whereas the one on the right is from calibration data. The 
vertical dashed line indicates periods of similar estimated exposure. 



iZip Monte Carlo iZip Calibration Data 
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Fig. 3 Longitudinal charge partition as measured in a CDMS iZIP detector. The partition is 
defined to be the hole signal minus the electron signal. The figure on the left shows a Monte 
Carlo response of the detector whereas the one on the right is from calibration data. The vertical 
dashed line indicates periods of similar estimated exposure. 



computed out to n = 27 and r = 81 constraining the solution to ^5 mm in the 
angular direction. 

The comparison between Monte Carlo and data do not show good agreement 
as seen in Figure [4] Shown in these figures is the response of the second to inner- 
most electrode to Am- 241 gamma rays. The Monte Carlo never shows a change 
in charge response despite being run for an order of magnitude longer than the 
charge degradation time scale seen in calibration data. 

Charges that stop at or near a detector face would result in a large amount of 
image charges in the metal film on the detector face. This results in a kernel in 
which charges stopped near a detector face have a small effect on electric field 
orientation. This would suggest that most of the degradation in the test device is 
due to charges trapping in the detector bulk and the lack of degradation in Monte 
Carlo is due to a charge trapping time that is too small for this detector. 



6 



-2 V bias Monte Carlo 



O 0.5 



0.5 1 1.5 2 2.5 3 

Integrated Charge Carriers [number] x 1 o 9 



-2 V Calibration Data 



Time [hours] 



Fig. 4 Response of ionization test device to Am- 241 gamma rays. The figure on the left shows 
a Monte Carlo response of the detector whereas the one on the right is from calibration data. The 
lines at high value, initially at ~200 in data are from gammas originating in the collimated source 
above the Q2 electrode. The vertical dashed line indicates periods of similar estimated exposure. 
The Monte Carlo results show an additional line with response ^0 for events generated above 
one of the other three electrodes. 



5 Conclusions and Future Work 



There is possible qualitative agreement between MC and calibration data in the 
iZIP radial charge partition and overall signal but poor agreement in longitudinal 
charge partition and the time scale for degradation in the iZIP and ionization test 
device. The Monte Carlo described in this paper did not take into account the 
increase in charge trapping rates at low fields and including this may improve the 
Monte Carlo results. Furthermore, the charge trapping times may be significantly 
higher in the ionization test device compared to what was included in Monte Carlo 
and this additional improvement in Monte Carlo may improve comparison. 
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